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Abstract The effects of global flow rotation and curvature on the subcritical transition to 
turbulence in shear flows are examined. The relevant time-scales of the problem 
are identified by a decomposition of the flow into a laminar and a deviation from 
laminar parts, which is performed for rotating plane Couette and Taylor-Couette 
flows. The usefulness and relevance of this procedure are discussed at the same 
time. By comparing the self-sustaining process time-scale to the time-scales 
previously identified, an interpretation is brought to light for the behavior of the 
transition Reynolds number with the rotation number and relative gap width in 
the whole neighborhood (in parameter space) of the non-rotating plane Couette 
flow covered by the available data. 

1. Introduction: 

In the last decade or so, a number of breakthroughs have been achieved in 
the understanding of the onset of turbulence in subcritical shear flows, such as 
the plane Couette flow and Poiseuille flow, both from an experimental point 
of view (e.g., Bottin et ai, 1998), and a numerical and theoretical one (e.g., 
Nagata, 1990; Clever and Busse, 1997; Hamilton et ai, 1995; Waleffe, 1997; 
Waleffe, 1998; Waleffe, 2003; Faisst and Eckardt, 2003). In this context, the 
present contribution has two main objectives: characterize the effects of global 
flow rotation and curvature in subcritical flows from the available data, and 
show that these characteristics can be understood at a semi-quantitative level 
from time-scale considerations. Understanding these questions is essential for 
geophysical and astrophysical applications, which is one of the motivations of 
this work. Data on rotating plane Couette flows and Taylor-Couette flows are 
used in this investigation. 

Section 1.2 is devoted to the identification of rotation and curvature char- 
acteristic quantities, and relating them to the gross dynamics of the flow. Not 
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surprisingly, the associated dimensionless numbers reduce to the shear-based 
Reynolds number, the rotation number, and the relative gap width (for Taylor- 
Couette flows); the novel point is that the relative gap width is interpreted in 
terms of a ratio of dynamically relevant time-scales. The experimental data 
are then reviewed in section 1.3, and the characterization of the available data 
in terms of the previously defined dimensionless numbers is performed in sec- 
tion 1.4. The last section briefly summarizes the most important results, and 
discusses their astrophysical implications. 

2. The physics of the advection term revisited: 

The main objective of this section is to pinpoint the relevant time-scales 
in globally subcritical, rotating and curved flows, and to relate them to the 
various contributions of the advection/acceleration term. This turns out to be 
essential to develop a semi-quantitative understanding of the available data on 
such flows. In practice, we consider only rotating plane Couette and Taylor- 
Couette flows. Incompressibility is assumed throughout. 

Equations of motions 

The relevant time-scales are well-known in rotating plane Couette flows, 
and follow immediately from the expression of the Navier-Stokes equation in 
the rotating frame, which reads 

du V7T 

— + u.Vu = 2ft x u + z/Au, (1) 

at p 

with obvious notations (in particular, the centrifugal term has been included in 
the pressure gradient). They are the shear 1 time-scale t s = IS' 1 ], the viscous 
one t v = cP/v (d is the gap in the experiment), and the rotation time-scale 
related to the Coriolis force tci = (2Q)~ l (Q, is the rotation velocity of the 
flow in an inertial frame), and relate to the advection term, the viscous term, 
and the Coriolis force, respectively. Correlatively, the flow is described by two 
dimensionless numbers, the Reynolds number 2 Re = t v /t s = \S\d 2 /u, and 
the rotation number Rq = sgn(S)t s /tn = 20,/ S. 

The situation is less straightforward for Taylor-Couette flows, where the 
dimensionless number usually associated to the flow geometry, 7] = ri/r Q (rj 
is the inner cylinder radius, r Q the outer one) does not obviously correspond to 
a ratio of time-scales of the flow. However, on closer inspection, it appears that 



'The convention adopted here is that the sign of S is chosen to be positive when the flow is cyclonic, i.e., 
when the contributions of shear and rotation to the flow vorticity have the same sign. With the usual choice 
of axes in plane Couette flows, this implies that S = —2S xy , where Sij is the usual deformation tensor. 
2 This definition differs from the usual one by a factor of 4; this convention is adopted here for consistency 
with the treatment of Taylor-Couette flows. 
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this situation arises from a partially incorrect assimilation of the shear time- 
scale to the advection term. Indeed, in the case of solid body rotation, the shear 
vanishes, while the advection term does not, due to the global curvature of the 
flow. One must therefore devise a way to isolate the global shear contribution 
to the advection term from other contributions. 

It turns out that one convenient way to operate such a distinction is to de- 
compose the total flow into its laminar part, and a (not necessarily small) devi- 
ation: 

u = u L + w. (2) 

Although the dynamical relevance of the laminar flow to the turbulent one is 
not a priori obvious, this procedure is suggested and justified by the following 
considerations: 

■ Inasmuch as this is feasible, a distinction between global shear, rotation 
and curvature cannot be operated by a tensorial decomposition of the 
advection term. For example, it is well-known that both a pure global 
shear and a global rotation, such as the ones present in rotating plane 
Couette flows, contribute to the vorticity tensor. In fact, a direct Taylor 
expansion of the deformation for small displacements shows that one 
needs to go at least to second order to distinguish the two contributions. 
Therefore, no tensor constructed from the flow velocity first derivatives 
will establish the required distinction, by construction. 

■ The global characteristics of the flow are the same for the laminar and 
turbulent solutions (geometry, global time-scales, nature of the bound- 
ary condition, etc). Therefore, one way to make them appear explicitly 
in the Navier-Stokes equation is precisely to make the proposed decom- 
position, as the laminar solution depends everywhere explicitly on these 
global characteristics. In particular, the laminar and turbulent flows share 
the same boundary conditions (velocity difference on the boundary, gap 
width, etc), so that the relative difference between the laminar and tur- 
bulent solution is of order unity. This means that the laminar solution is 
a convenient measure of the turbulent one, although their detailed mech- 
anisms and characteristics are of course essentially different. For exam- 
ple, it turns out that the transition Reynolds number is highly sensitive 
to various global and/or qualitative characteristics of the laminar flow, 
such as time-scales, or "distance" in parameter space to the linear stabil- 
ity limits (see section 1.3). 

It is useful to point out where this decomposition leads to for the rotating 
plane Couette flow. This is most naturally done in the rotating frame, so that 
Eq. (1) becomes 
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<9w <9w V(57r 

— + w.Vw = S ■ y— + (20 + - 2Clw x e v h z^Aw, (3) 

at ox p 

where the pressure gradient balancing the laminar flow Coriolis force has been 
subtracted out to form the effective generalized pressure 5ir. Note that on 
the walls, which specify the global characteristics of the flow, the boundary 
condition becomes w = 0, quite a featureless constraint: the effect of these 
boundaries is now explicitly included in the Navier-Stokes equation through 
the dynamical linear forcing terms on the right-hand side. The real usefulness 
of this change of point of view comes out when considering Taylor-Couette 
flows, as we shall now argue. 

In this flow, the laminar solution takes the form U£,(r) = rCl(r)eg. Note that 
the rotating plane Couette flow can be viewed as a limit of small relative gap 
width d/f — > (f is some characteristic radius of the flow), at constant shear, 
and constant rotation. Then, the effect of the global curvature and rotation will 
be more easily distinguished from one another if one chooses a formulation of 
the Navier-Stokes equation which makes the difference with Eq. (3) explicit. 
To this effect, one must define a characteristic rotation velocity O, and a char- 
acteristic shear 3 S of the laminar flow. A convenient way to do this is to choose 
a characteristic radius f, and impose Q = 0(f), S = S(r); the choice of f does 
not need to be further specified for the time being (this point is discussed in the 
next subsection). Defining 5£l = 0(r) — and 5S = S(r) — S, the decompo- 
sition of the Navier-Stokes equation of the Taylor-Couette flow leads to (r, 4>, z 
is the coordinate system in the rotating frame) 



— + w.Vw = - 5Q,d' 6 w - (2Q + SWe,* + 2Qw (t) e r h Mw 

at v p 

- {28VL + dS)w r e,p + 2(Klu^e r , (4) 

where one has used <9^w = d^-w+5^}w r e ( p+2S^}w ( f ) e r , with <9^w = (d ( j > w r )e r + 
(dtpw^ecf, + (d(j,w z )e z ; this definition is introduced so that the contributions 
of order 1 jf of the derivatives in the linear terms are separated from the con- 
tributions of order 1/d. This equation is similar 4 to Eq. (3), except for the 
last two terms, which consequently are connected to the global flow curvature. 
Note that, although the definition of 5£l and 5S depends on the choice of f, the 
overall variation of these quantities throughout the flow, A(<5$7) and A(5S), 



3 The shear of the laminar flow is defined as S = rdU jdr = 25 r ^ in order to maintain the sign convention 
adopted for rotating plane Couette flows for cyclonic and anticyclonic rotation. 

4 In the identification of the two equations, note that r < > y and <f> < > — x. Also, 8Q ~ S(r — f)/f, 

an approximation which holds to ~ 10% for the range of j) explored in the available experiments, a feature 
needed in the comparison of Eqs. (3) and (4). 
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does not. In fact, |A(<5fi)| ~ |A(<5S)| ~ \S\d/f, as can be checked from 
the laminar flow profile f2(r) = A + B/r 2 . In the process, four time-scales 
of the incompressible Taylor-Couette flow have been identified: they are the 
shear time-scale t s = IS" 1 ], the rotation time-scale tn = (20) _1 , a curvature- 
related time-scale that one can define as tc = |<S _1 |f/d, and, of course, the 
viscous time-scale t v = d? /v. This also shows that the dimensionless geomet- 
ric number rj can be related to the ratio of the shear and curvature time-scales. 
Taking the limit d/f — ► (rj — > 1), one recovers the rotating plane Couette 
relation Eq. (3). 

Characteristic quantities, dimensionless numbers, and the 
curvature and rotation concepts: 

Taylor-Couette flows possess three dimensionless control parameters, which 
are usually chosen as the Reynolds numbers associated to the inner and outer 
cylinder rotation velocities, Ri and R Q , and rj, the ratio of their two radii. The 
preceding discussion suggests that one should use instead ratios of time-scales, 
which have a more direct dynamical meaning. This defines the shear-based 
Reynolds number 5 Re = t v /t s , the rotation number Rq = sgn(S)t s /tn, and 
a "curvature" number Rc = t s /tc- Once a choice of f is operated (see below), 
the procedure is completely specified. This three-dimensional parameter space, 
in which notable curves and surfaces are drawn, is represented in figure 1. 

However, the physical meaning of this procedure is less straightforward than 
one would like, and this is related to an obvious weakness of Eq. (4): the 
distinction between the rotation and curvature terms is not absolute when the 
related time-scales are both dynamically significant. Indeed, any change of 
definition of Cl and S results in a correlative change of SQ and 5S. Never- 
theless, the physical meaning is to a large extent unambiguous in at least two 
different contexts: 

■ If one changes the rotation velocity of the inner and outer cylinders by 
the same quantity, this will change f2 by the same amount (indepen- 
dently of the choice of f), while leaving all other quantities (S, 5Q, 6S) 
unchanged. Such changes are obviously an effect of changes in the flow 
rotation. 

■ On the other hand, when changes in the flow are operated while main- 
taining tn < tc, the physical meaning of the distinction between the 
characteristic quantities (fj, S) and the deviations from these (SCI, 5S) 



5 We use the same Reynolds number definition for rotating plane Couette and Taylor-Couette flows, based 
on the total gap width and total velocity difference. Consequently, the quoted Reynolds numbers for plane 
Couette flows differ from the ones in the literature by a factor of 4. 
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is blurred. This is the case in particular when the cylinders are counter- 
rotating, or when one cylinder is at rest, for any choice of f. In such a 
context, changes in both parameters (the rotation number and the cur- 
vature number) describe the effect of a change in the flow curvature, as 
they are both proportional to d/f, and vanish if the limit of vanishing 
global curvature is taken while enforcing the tn < tc relation. 

This clearly shows that rotation and curvature are not interchangeable con- 
cepts, although they have a non-negligible overlap. In this context, the denomi- 
nations "rotation number" and "curvature number" are somewhat conventional 
and partially misleading, even if justified to some extent by the preceding con- 
siderations. In the (Rq, Rc) plane, changes along lines of constant Rq cor- 
respond to changes of rotation, but there are infinitely many paths involving 
changes of both parameters and corresponding to changes of curvature from 
a physical point of view. Furthermore, most paths do not lead to any clear- 
cut distinction between curvature and rotation changes in the flow. Of course, 
once a reference curvature path is chosen (e.g., the path with the inner cylin- 
der at rest, corresponding to the inviscid linear stability limit on the cyclonic 
side, shown as curve (4) on figure 1), every point in the (Rq, Rc) plane can be 
connected to the non-rotating plane Couette flow (the origin in the plane), first 
through a change of curvature along the chosen curvature path until the desired 
curvature number is reached, and then through a change of rotation at constant 
Rc- However, this distinction is only relative. 

Obviously, this situation is intrinsic, as one cannot curve a straight flow, 
without at the same time making it rotate. The procedure outlined here nev- 
ertheless leads to the definition of well-defined parameters, which have a dy- 
namical interpretation. It is their physical meaning in terms of rotation and 
curvature which is partially ambiguous. Furthermore, these parameters turn 
out to be useful to understand basic features in the data on the subcritical tran- 
sition to turbulence, as discussed in the next section. 

The remaining point to be addressed relates to the choice of f. The preced- 
ing discussion makes it clear that this choice is not unique. The definition we 
have adopted here is f = (r^ro) 1 / 2 , as suggested in Dubrulle et al, 2004 (r« 
and r are the inner and outer radius, respectively). This choice is partially mo- 
tivated by the compactness of the resulting expressions for the dimensionless 
numbers introduced above: 



Re = 



2f f\tt -tti\d 
r + ri v 



1 + 7? 



2 



\t]R q - Ri 



1 1 > 



(5) 



Rn = 



rjVLi + r Q d 
T Ti Vl 



(1-7?) 



Ri + Ro 

i]R - Ri ' 



(6) 
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Note that, for the range of values of 77 explored in the available experiment, 
Re ~ f|Afi|d/i/, and ifo ~ (2fi/f)(d/Afi) (within a few percents). 

3. Subcritical transition in rotating plane Couette and 
Taylor-Couette flows: 

Considering a laminar flow of given dimensionless numbers (rotation, cur- 
vature, and Reynolds), two different things can happen when increasing the 
Reynolds number: either the flow will undergo a linear instability first (super- 
critical transition), or it will undergo a laminar- turbulent transition first (glob- 
ally subcritical transition). The second option may happen whether the flow is 
linearly unstable or not. 

The Reynolds number characterizing subcritical transition in a system is not 
a quantity that can be measured with absolute precision, as it depends to some 
extent on the experimental protocol used in its determination. For example, 
the laminar-turbulent transition Reynolds number is genetically larger than the 
turbulent-laminar one. Furthermore, the flow is intermittent over a range of 
Reynolds numbers in the vicinity of this transition. This leads to some differ- 
ences in the determined Reynolds transition values, even when the same data 
are used by different authors; however, the dispersion of the data in a given au- 
thor's choice is much smaller. Overall, the resulting range of values (at given 
dimensionless numbers) is uncertain within less than a factor of ~ 2; we shall 
ignore this problem here, as we are only interested in characterizing qualitative 
trends and orders of magnitude. 

With this convention, both supercritical (Re = R c ) and globally subcritical 
(Re = R g ) transitions are characterized by surfaces in the three dimensional 
space (Rq,Rc, Re). Only particular lines on these surfaces have been probed 
by the available experiments. Obviously, the supercritical and subcritical sur- 
faces meet somewhere in this space, so that one needs to characterize both 
surfaces. 

It turns out, from a practical point of view, that the supercritical surface 
(manifold (1) in figure 1) is sufficiently well captured by the analytic approx- 
imation derived by Esser and Grossmann, 1996, for Taylor-Couette flows, al- 
though it relates only to axisymmetric perturbations (in other words, the non- 
axisymmetric perturbations seem to play little role in the definition of the su- 
percritical transition). The rotating plane Couette flow is included in the limit 
r\ — > 1. For the relatively high Reynolds numbers of interest for subcritical 
transitions to turbulence (> 1000), the dependence of the supercritical surface 
on the Reynolds number is very steep, and the surface is well-approximated 
by the inviscid linear stability limit (curves (2) of figure 1). This explains that 
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Figure 1. Stability boundaries in the (Re, Rn, Re) parameter space (77 is used as a measure 
of Rc on this graph). Manifold (1): linear stability threshold. Dark curves (2): linear stability 
threshold in the inviscid limit. Manifold (3): extrapolation of the inviscid criteria throughout the 
full Reynolds space (partially shown for readability). Curve (4): globally subcritical threshold 
obtained with inner cylinder at rest (see text). Other curves: globally subcritical thresholds 
obtained at fixed value of 77 [(5): 77 = 0.7; (6): 77 = 1] (see text). 



the limit of the subcritical regime in the (Rq, Rc) plane is well-approximated 
by the inviscid limit. Linear instability follows somewhere in the fluid in this 
limit if 

2Q(r) 
S{r) 

at this location (this is equivalent to the constraint put by the Rayleigh discrim- 
inant). Asking that the fluid is everywhere stable with respect to this criterion 
translates into Rq < —1, or Rq > (1 — rf)/r]. 

The data on subcritical transition discussed here are those of Wendt (1933), 
Taylor (1936), Tillmark and Alfredsson (1996), and Richard (2001); they are 
represented on figure 2. One could also include data on counter-rotating cylin- 
ders (Andereck, Liu and Swinney, 1986; Prigent et ai, 2003), but these occupy 
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Figure 2. Data on subcritical transition. Upper left: data of Wendt (1933) and Taylor (1936), 
for cyclonic rotation, with the inner cylinder at rest, for varying r\; this shows the dependence 
on the curvature number in cyclonic flows. Upper right: data of Tillmark and Alfredsson (1996) 
for rotating plane Couette flows {r\ = 1). Lower left and lower right: data of Richard (2001), 
for anticyclonic and cyclonic rotation, respectively (77 = 0.7). On the last three graphs, the 
curvature number is held fixed, so that the data show the dependence of the transition Reynolds 
number on the rotation number. The solid lines represent best quadratic or linear fits to the data. 
— (1 — v)/ 7 ! f° r cyclonic flows, and Rq = —1 for anticyclonic ones. 



only a small area in the (Rn, Rq) plane, and bring little constraint on the trend 
of the transition Reynolds number with the rotation and curvature numbers (see 
Longaretti and Dauchot, 2004 for a discussion o couter-rotating data). 

The data are well parameterized by the following approximate formula (the 
+ and — sign refer to cyclonic and anticyclonic flows, respectively) 
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R±(Ra,Rc) = R$ c + a^WlRn ~ Rn^l + b ± R 2 c , (9) 
with R+ c ~ 1400, R PC ~ 1100 ~ fl£ c , 21000 < a* < 61000, 2 x 10 5 < 
6+ < 6 x 10 5 , < 6+, and where R™' + = (1 - r/)/r/, and ~ = -1. 
This is discussed in detail in Longaretti and Dauchot, 2004, and Dubrulle et 
al, 2004. 

The most notable characteristics of this dependence are the following: 

■ The linear dependence on the rotation number, and quadratic one on the 
curvature number. 

■ The rather steep dependence with both numbers (a ± and b + are large 
numbers). 

■ The apparent symmetry between cyclonic and anticyclonic rotation num- 
ber dependence, and the absence of dependence on the curvature on the 
anticyclonic side (b~ ~ 0). 

The next section discusses these features; in particular, the strength of this 
dependence is explained in order of magnitude on the basis of time-scale con- 
siderations. 

4. Data interpretation: 

The linear and quadratic dependence just pointed out can be viewed as the 
result of a Taylor expansion. For cyclonic rotation, this expansion is performed 
around the non-rotating plane Couette flow, first in terms of Rc along curve 
(4) on figure 1, and then away from it, at constant Rc, in terms of Rq. For 
anticyclonic rotation, the expansion depends only on Rq, at least for r\ > 0.7, 
and is performed from the marginal stability state, at Rq ~ — 1. 

Linear dependence 

The linear dependence with Rq is neater on the cyclonic data than on the an- 
ticyclonic one, but the data extend less far on the cyclonic side, unfortunately. 
The increased dispersion on the anticyclonic side has several reasons: 

■ Both cylinders need to be rotating at much higher speed than on the 
cyclonic side to reach the subcritical turbulent transition. This automat- 
ically reduces the precision of the measurements. 

■ The quantity Rq amplifies the uncertainty due to small errors in the de- 
termination of the cylinder angular velocities at transition. 

■ There is an important intrinsic difference between the cyclonic and an- 
ticyclonic marginal stability limits. On the cyclonic side, instability be- 
gins at a single location (the inner radius), whereas on the anticyclonic 
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side, marginal stability applies throughout the fluid (this follows because 
the constant angular momentum solution is a laminar solution of the 
Taylor-Couette flow). Equivalently 250. (r) + SS(r) = at the anticy- 
clonic marginal stability limit. Therefore, the fluid is much more sensi- 
tive to a potential linear instability. In particular, the unavoidable Eck- 
mann circulation can much more easily make the criterion Eq. (8) satis- 
fied somewhere in the flow on the anticyclonic side than on the cyclonic 
one. We believe that this feature most likely explains why Richard's 
data show only a weak dependence of the transition Reynolds number 
on the rotation number out to 1 1 + Rq | ~ 0.35, and then a sharp increase 
to reach back the linearly varying regime. The related data points are 
shown as thin dots in the lower left quadrant of Fig. 2, and not used in 
the linear fit. 

Note that the mutual cancellation of a part of the curvature terms in Eq. (4) at 
the anticyclonic marginal linear stability limit just pointed out also provides a 
natural explanation for the apparent absence of dependence of the anticyclonic 
data on the curvature number. 

Quadratic dependence 

The quadratic dependence of the data on the gap width when keeping the 
inner cylinder at rest has already been pointed out by Zeldovich, 1981, and 
Richard and Zahn, 1999. At least three explanations of this behavior have been 
put forward in the literature (Zeldovich, 1981; Dubrulle, 1993; Longaretti, 
2002). These will be commented below. 

Let us first note that, in plane Couette flows, transition occurs at a constant 
Reynolds number Re = Sd 2 jv. Therefore, at a given shear S, the scale d is the 
characteristic scale of turbulence in such a system. Conversely, in the quadratic 
asymptotic regime, [Re oc (d/f) 2 ], transition occurs at constant Re* = Sf 2 jv 
(= b + ), a point already made by Richard and Zahn, 1999. Consequently, the 
characteristic radius f instead of the gap d characterizes the transition at a given 
shear. This unambiguously points out curvature and not rotation as the source 
of this behavior, consistently with the discussion of section 2. 

Dubrulle, 1993, explains the quadratic behavior by considering the growth 
of finite amplitude local defects in the laminar profile. However, only WKB 
modes of instability created by the defects are considered, for which the scale 
f cannot play any role. This is why we believe that this analysis cannot capture 
the transition mechanism. 

Zeldovich, 1981, phenomenological explanation is based on the following 
two ideas: 

■ The transition Reynolds number may depend on the single time-scale 
ratio Ty(r) = An 2 (r) / S 2 (r) at some appropriately chosen radius in the 
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flow: R g = f(Ty); K (r) = [20(20 + S")] 1 / 2 is the epicyclic frequency, 
i.e., the frequency of oscillation of the whole fluid, under the combined 
action of the shear and the Coriolis force. 

■ A "split-regime" of instability may occur, in which the inner portion of 
the flow undergoes a transition to turbulence at lower Reynolds number 
than the whole flow, once a large enough relative gap width is reached. 

Considering the relative ambiguity in the distinction between the rotation 
and curvature time-scale discussed in section 1.2, one may indeed ask whether 
a single time-scale would be sufficient to understand the data. This would 
imply that Taylor-Couette flows possess a hidden redundancy, and could be 
described by two appropriately chosen parameters, and not three. However, 
the extended set of data used here does not support this idea, although a larger 
body of experimental results is probably needed to ascertain this result. 

Longaretti, 2002 has developed an alternate phenomenology of subcriti- 
cal shear flow turbulence. It relies on the one hand on a turbulent viscosity 
description, in which the characteristic length is identified to the top of the 
Kolmogorov cascade; on the other hand it makes use of the constraint that 
a subcritical shear flow is out of thermodynamic equilibrium and tries to re- 
store equilibrium by transporting momentum across the shear, and, to do so, 
chooses the most efficient of the two means at its disposal (laminar or turbulent 
transport). This provides scaling laws for the characteristic length and velocity 
of turbulent eddies, by relating them to the Reynolds number of transition to 
turbulence. By further noting that for large enough gap widths, the turbulent 
eddies must unavoidably scale with the radius and not the gap, the quadratic 
regime is recovered. 

Orders of magnitude 

The last point we wish to address is the origin of the large values the co- 
efficients a and b which appear in the characterization of the data performed 
in Eq. (9). We mostly consider cyclonic flows, for which the physics is best 
understood 

As mentioned in the introduction, an important breakthrough in the under- 
standing of subcritical turbulence in non-rotating plane Couette flow comes 
from the work of Waleffe, 1997, who analyzed by the means of quasi-linear 
theory a turbulent self-regeneration process previously observed in the numer- 
ical simulations of Hamilton et al, 1995. These last authors have tracked down 
turbulence to the smallest unit where it is self-sustained by reducing appropri- 
ately the simulation box size and the Reynolds number. A very important 
feature of the identified self-sustaining process is that it has a rather long time- 
scale compared to the shear: 
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t+^ioos- 1 . (10) 

and that the scales involved in the self-regeneration mechanism are compara- 
ble to the flow width. This time-scale is the shortest of all the mechanisms 
found in the systematic reduction of the flow, and thus corresponds to the most 
robust one, which involves two streamwise rolls in the spanwise direction. 
These streamwise rolls, first observed by Dauchot and Daviaud, 1995, typi- 
cally scale on the gap width. Accordingly, it is very likely that such a long 
time-scale is a generic feature in non rotating plane Couette flows, because of 
the large (Rp C ~ 1500) Reynolds number of transition to turbulence which 
are always observed in these systems. Such a large Reynolds number con- 
strains the viscous diffusion time at the scales involved in the self -regeneration 
mechanism. Typically for a length-scale d/A, the equality of the viscous time 
t v = {d/4) 2 /v with i+ would indeed leads to R g = 1600. Such a scale 
(d/4) is characteristic of the thickness of the streaks, apparently the smallest 
characteristic scale of the process. 

From a physical point of view, one expects that the contributions of either 
rotation or curvature become comparable to Rp C in Eq. (9) when the rotation 
or curvature time-scales decrease to become comparable to t+ . Correlatively, 
the fact that this happens when the rotation and curvature numbers are ~ 1/20 
comes as no coincidence. Indeed, the timescale associated with the Coriolis 
term is ta ~ (2S7) -1 , while the timescale associated with the curvature terms 
is tc — (<5f2) -1 ~ Sf/d. They become comparable to tf sp as defined in 
Eq. (10) when Rq or Rc exceeds 10~ 2 , which is remarkably close to the actual 
value of 1/20, considering the qualitative nature of the argument. This physical 
constraint is also what primarily determines the magnitude of a + and b + (once 
the form of the dependence on Rq and R c is known). Indeed, in rotating plane 
Couette flows, requiring o, + Rq > R e t>c wnen Rn ^ 1/20, implies that 
a + ~ 10 4 . Similarly, in Taylor-Couette flows, requiring that b + R 2 > Re~p C 
when R c > 1/20 leads to b+ ~ 10 5 . 

One can see that the rather large values of the transition Reynolds num- 
ber, as well as of the coefficients characterizing the effect of rotation (a + ) and 
curvature (6 + ) can be ascribed to a single origin: the rather large ratio of the 
turbulence self-regeneration process to the shear time scale. Still, we do not 
infer that the self-sustaining mechanism proposed by Waleffe is valid in the 
presence of rotation. We just use the fact that in order to modify or even sup- 
press this mechanism, the rotation or curvature effects must have timescales of 
the same order. On the contrary, the above analysis clearly indicates that a bet- 
ter understanding of this process in the framework of curved shear flows, and 
identifying it in the presence of rotation, is of primary importance for future 
progress. 
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The self-sustaining process has not yet been identified at the anticyclonic 
marginal stability limit, and its nature is not known. However, one can reverse 
the reasoning expressed right above to reach the conclusion that its character- 
istic time-scale is also ~ lOOS -1 . Testing this conjecture would bring support 
to the framework developed here to analyze and understand the data. 

5. Conclusions and implications: 

We feel that this work brings to light a few important points, which we 
believe to be of potentially more general applicability than what was done 
here: 

■ It is both meaningful and useful to decompose the Navier-Stokes equa- 
tion into a laminar part of the flow, and a deviation from the laminar flow. 
This helps identifying the relevant time-scales in the flow, and isolating 
which "portion" of the advection term is directly related to the physics 
described by the Reynolds number. The same procedure is also useful to 
relate various flows to one another. 

■ The general trends and features of the transition Reynolds number data 
can be understood with the help of the previously identified time-scales, 
and over a reasonably extended fraction of the parameter space, once the 
turbulence self-sustaining process time-scale is identified at one point in 
the parameter space. The procedure in turns constrains to some extent 
the time-scale of the self-sustaining process in the domain where data 
are available, a point we have only briefly touched upon in the preceding 
section. 

In the process, we have tried to elucidate a little more the relation between 
the flow global rotation and its global curvature. We found that rotation and 
curvature effects most probably cannot be accounted for by a single time-scale, 
as suggested by Zeldovich, 1981, but that two time-scales are required, in ac- 
cordance with the fact that Taylor-Couette flows are characterized by two di- 
mensionless numbers besides the Reynolds number. 

The data described here, and the analysis we have developed, have some 
bearing to a related astrophysical problem, namely, the existence of subcritical 
turbulence in keplerian accretion disks, a question which has been the object 
of an important debate in the astrophysical community over the three or four 
past decades. Such disks are observed in relation to the formation of young 
stars; they are also believed to be present in a variety of other astrophysical 
objects, such as active galactic nuclei, microquasars, and so on. A large scale 
keplerian profile (U(r) oc r~ 3 / 2 ) follows if the disk is cold enough. The profile 
can nevertheless stochastically deviate from keplerian on scales comparable to 
the disk scale height h <C r. Young (proto-)star disks are also probably not 
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ionized enough for MHD processes to be relevant over a significant fraction of 
their extent. 

The microscopic transport in these disks is known to be many orders of 
magnitude smaller than the one inferred from the observations, so that these 
disks are believed to be turbulent. Hydrodynamic Keplerian disks are linearly 
stable in their most simple flavors. It was believed until the mid-90's that they 
were nevertheless hydrodynamically turbulent, on the basis of the experimental 
evidence of subcritical turbulence in non-rotating Couette flows. 

This belief was challenged by the numerical simulations of Balbus et at, 
1996, and Hawley et ai, 1999. These authors have studied whether keple- 
rian disks are locally turbulent, by reducing the Navier-Stokes equation in the 
disk to Eq. (3), with "shearing-sheet" boundary conditions (they ignore the 
disk vertical stratification), thus asking the question in a framework which is 
extremely close to the one studied here. They found that a dynamically signif- 
icant and stabilizing Coriolis force prevents the appearance of turbulence for 
keplerian-like flows (Rq = —4/3), up to to the highest resolution achieved in 
the simulation (256 3 with finite-difference codes; see the referenced papers for 
detail). 

These simulations have had a very large impact in the astrophysics commu- 
nity, where the now most largely spread opinion is that a linear instability is 
needed for turbulence to show up, a somewhat excessive position. However, 
the experimental results of Richard, 2001, seem to imply that keplerian-like 
flows should be turbulent at very modest (in astrophysical standards) Reynolds 
numbers. The problem is currently reinvestigated from a numerical point of 
view (Lesur and Longaretti, in preparation). Preliminary results indicate that 
the Coriolis force does not prevent the existence of turbulence in rotating plane 
Couette flows, but of course alters the turbulence properties, as one would 
expect. However, the relevance of subcritical hydrodynamic turbulence to ac- 
cretion disk transport remains to be more precisely investigated. 

References 

Andereck, CD., Liu, C.C., and Swinney, H.L. (1986). Flow regimes in a circular Couette system 
with independently rotating cylinders. J. Fluid Mech., 164: 155-183. 

Balbus, S.A., Hawley, J.F., and Stone, J.M. (1996). Nonlinear stability, hydrodynamical turbu- 
lence, and transport in disks, Astrophys. J., 467: 76-86. 

Bottin, S., Dauchot, O., Daviaud, F., and Manneville P. (1998). Experimental evidence of stream- 
wise vortices as finite amplitude solutions in transitional plane Couette flow. Phys. Fluids, 
10, Issue 10: 2597-2607. 

Clever, R.M., and Busse, FH. (1997). Tertiary and quaternary solutions for plane Couette flow. 
J. Fluid Mech., 344: 137-153. 

Dauchot, O., and Daviaud, F. (1995). Finite amplitude perturbation and spot growth mechanism 
in plane Couette flow, Phys. of Fluids 7: 335-343. 



16 



Dubrulle, B. (1993). Differential rotation as a source of angular momentum transfer in the Solar 

Nebula, Icarus, 106: 59-76. 
Dubrulle, B., Dauchot, O., Daviaud, R, Longaretti, P.-Y, Richard, D., and Zahn, J.-P. Stability 

and turbulent transport in rotating shear flows: prescription from analysis of cylindrical and 

plance Couette flow data. Submitted to Phys. of Fluids. 
Esser, A., and Grossmann, S. (1996). Analytic expression for Taylor-Couette stability boundary, 

Phys. Fluids, 8: 1814-1819. 
Faisst, H., and Eckardt, B. (2003). Travelling waves in pipe flow. Phys. Rev. Lett., 91: 224502. 
Hamilton, J.H., Kim, J., and Waleffe, F. (1995). Regeneration mechanisms of near-wall turbu- 
lence structures. / Fluid Mech., 287: 317-348. 
Hawley, J.F., Balbus, S.A., and Winters, W.F. (1999). Local hydrodynamic stability of accretion 

disks, Astrophys. J., 518: 394-404. 
Longaretti, P.-Y. (2002). On the phenomenology of hydrodynamic shear turbulence, Astrophys. 

J., 576, 587-598. 

Longaretti, P.-Y., and Dauchot, O. Global rotation-curvature time-scales, and the subcritical 
transition to turbulence in shear flows. Submitted to Phys. of Fluids. 

Nagata, M. (1990). Three-dimensional finite-amplitude solutions in plane Couette flow: bifur- 
cation from infinity. J. Fluid Mech., 217: 519-527. 

Prigent, A., Gregoire, G., Chate, H., and Dauchot, O. (2003). Long wavelength modulation of 
turbulent shear flows. Physica D, 174: 100-1 13. 

Richard, D. (2001). lnstabilites hydrodynamiques dans les ecoulements en rotation differentielle. 
PhD thesis. Universite de Paris VII. 

Richard, D., and Zahn, J.-P. (1999). Turbulence in differentially rotating flows, Astron. Astro- 
phys., 347: 734-738. 

Taylor, G.I. (1936) Fluid friction between rotating cylinders, Proc. Roy. Soc. London, A 157: 
546-564. 

Tillmark, N., and Alfredsson, PH. (1996). Experiments on rotating plane Couette flow. In Ad- 
vances in Turbulence VI. (eds, Gavrilakis et ah), p. 391-394; Kluwer. 

Waleffe, F. (1997). On a self-sustaining process in shear flows. Phys. Fluids, 9: 883-900. 

Waleffe, F. (1998). Three-dimensional coherent states in plane shear flows. Phys. Rev. Lett., 
81:4140-4143. 

Waleffe, F. (2003). Homotopy of exact coherent structures in plane shear flows. Phys. Fluids, 
15:1517-1534. 

Wendt, G. (1993). Turbulente Stremung Zwischen Zwei Rotierenden Konaxialen Zylindern, 
Ing.Arch.,4: 577-595. 

Zeldovich, YB. (1981). On the friction of fluids between rotating cylinders, Proc. Roy. Soc. 
London A, 374: 299-312. 



